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Abstract. The coming years will see routine use of solar data of unprecedented spatial and 
spectral resolution, time cadence, and completeness in the wavelength domain. To capitalize on 
the soon to be available radio facilities such as the expanded OVSA, SSRT and FASR, and the 
challenges they present in the visualization and synthesis of the mult i- frequency datasets, we 
propose that realistic, sophisticated 3D active region and flare modeling is timely now and will 
be a forefront of coronal studies over the coming years. Here we summarize our 3D modeling 
efforts, aimed at forward fitting of imaging spectroscopy data, and describe currently available 
3D modeling tools. We also discuss plans for future generalization of our modeling tools. 
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1. Introduction 

Solar activity, although energetically driven by subphotospheric processes, depends 
critically on coronal magnetism, which, broadly speaking, includes magnetic field gener- 
ation, evolution, and transformation into kinetic, thermal, and nonthermal energies in the 
corona. Reliable tools for doing direct diagnostics have been lacking, although the situa- 
tion is currently changing. Indeed, new space- and ground-based solar optical telescopes 
are already capable of precise measurements of the photospheric magnetic field with sub- 
arcsecond angular resolution and high temporal resolution. When combined with modern 
extrapolation algorithms, these data offer important clues on the coronal magnetic field 
structure and evolution. However, given the finite angular resolution, sensitivity, observa- 
tional errors, and even theoretical limitations, such extrapolations are not unique, so the 
extrapolations require independent verification. An opportunity for quantitative verifi- 
cation will be available when the new generation of high-resolution solar-dedicated radio 
instruments (expanded OVSA, SSRT, and FASR) become operational. Microwave radi- 
ation is produced by the gyrosynchrotron (GS) mechanism as accelerated fast electrons 
gyrate in the coronal magnetic field. As has been recently proven using simulated mi- 
crowave data, the coronal magnetic field can indeed in principle be reliably recovered at 
the flare dynamic time scales from the radio data, along with the key parameters of the 
thermal plasma and accelerated electrons (Fleishman et al. 2009). The ability to detect 
the magnetic field and its changes on dynamic time scales is a critically needed element 
to uncover the fundamental physics driving solar flares, eruptions, and activity. 



2. Modeling Methods and Tools 

Direct Modeling. 3D models of the solar flares are not yet numerous. Those available 
models (Preka-Papadema & Alissandrakis 1992, Kucera et al. 1993, Bastian et al. 1998, 
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Figure 1. Left: interface of the model geometry and non-uniform parameter input — model mag- 
netic loop (left window); nonuniform distribution of fast electrons given in a color code (middle 
window); and numeric input and graph display of the parameters involved. Right: selection of the 
fast electron energy distribution from a number of pre-defined analytical functions and parameters 
of this distribution. 



Simoes & Costa 2006, 2010, Tzatzakis et al. 2008, Fleishman et al. 2009) are built on an 
idealized (e.g., dipole) magnetic loop, rather than a realistic magnetic geometry. 

Our currently available modeling tool, GS Simulator, is also built based on an an- 
alytical (dipole) magnetic field model. A flaring loop model (i.e., user specified dipole 
fiux tube. Figure [U left) is produced with a newly developed interactive IDL widget 
application intended to provide a fiexible tool that allows the user to generate spatially 
resolved GS spectra. To do so, the user populates the loop with thermal plasma by se- 
lecting plasma temperature and density in the adjustable parameter list. Figure [H right, 
and specifies a fast electron population by selecting one of a few pre-defined distributions 
of the electrons over energy and pitch-angle and choosing numeric parameters for these 
distribution (Figure [TJ right). The spatial distribution of the electrons is also specified 
along the loop, as shown in the color coded image (Figure [TJ left). After selecting the 
angle from which to view the loop with the mouse, the tool then calculates the physi- 
cal parameters along each line of sight needed to solve the radiation transfer equation, 
which is performed by external callable computing blocks. The default codes generat- 
ing the GS (and free-free) emission based on the input geometrical line-of-sight model 
data were written in FORTRAN and C++ based on fast GS codes newly developed by 
Fleishman fc Kuznetsov (2010)[ and compiled as a DLL (or SO in case of Linux) callable 
by IDL. 

The object-based architecture of this application provides the user with full 3D inter- 
action with a predefined, but adjustable, magnetic loop geometry, as well as with any 
user defined analytical geometrical model that would inherit the basic properties of the 
generic "gs.model" IDL object defined in this package. To allow for more realistic 3D 
fiare modeling, in place of drawing a loop "by hand", the tool must be further devel- 
oped to include (i) a numerical magnetic field structure, such as would be obtained from 
a photospheric extrapolation (Figure [2]) or a full MHD model based on vector photo- 
spheric measurements of the magnetic and velocity fields, along with the thermal plasma 
distribution, (ii) a realistic electron acceleration and transport model, and (iii) the abil- 
ity to quickly compute emission in various wavelength regimes. The tool will enable the 
user to identify, from the model or through comparison with observations, the subset of 
field lines involved in fiaring, consider fast electron transport in this realistic magnetic 
structure, calculate radiation from this evolving volume, and so simulate a solar fiare. 

Forward Fitting. The coronal magnetic field is a key parameter controlling most 
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Figure 2. Left: Two flux tube models. They can be defined independently, and may represent 
interacting loops. Additional controls are being prepared for populating the loops with particles and 
ambient medium. Even dynamics of the particles is planned. Right: The loops now viewed from a 
specified location on the solar disk (22 E 21 S in this example). The tabs labeled Flux Tube 3 and 
Flux Tube 4 represent the two loops. 



solar flaring activity, particle acceleration and transport. It has been understood, and 
often proposed, that the coronal magnetic field along with fast electron distribution can 
in principle be evaluated from the microwave GS radiation, which is indeed sensitive to 
the instantaneous magnetic field strength and orientation relative to the line of sight and 
to the fast electron spectrum. 

Diagnostics, understood as the determination of physical parameters of a system under 
study from arrays of observed parameters, is a key outstanding problem in Solar Physics. 
In some (basically linear) cases regularized true inversions can work well (e.g., Kontar et 
al. 2004). In most of the cases, however, such true inversions fail because of the highly 
nonlinear nature of physical systems. In such cases, the forward fitting, i.e., finding a 
number of free parameters of a physically motivated model of the system from fitting the 
model to observations, can often be used successfully in place of true inversions. 

Anticipating a large breakthrough in the radio imaging spectroscopy observations, 
which will become possible soon due to the next generation of radio instruments, we 
have developed a practical forward fitting method, based on the SIMPLEX algorithm 
with shaking, that allows reliable derivation of the magnetic field and other parameters 
along a solar flaring loop using microwave imaging spectroscopy of GS emission, which 
is calculated with newly developed fast GS codes (Fleishman & Kuznetsov 2010). We 
illustrate the method using a model loop with spatially varying magnetic field, filled 
with uniform ambient density and an evenly distributed fast electron population with an 
isotropic, power-law energy distribution (Fleishman et al. 2009). 

From the flare radio model (Figure [3l left) described above we have a sequence of 
spatially resolved microwave spectra (in a general case, both intensity and polarization 
data, one spectrum per pixel). Then, we fit the data to a model microwave spectrum 
pixel by pixel (Figure [3l middle) to derive physical parameters of the source (e.g., the 
magnetic field, Figure[3l right). Although the exact GS formulae are very computationally 
expensive, much faster codes giving the same accuracy have recently been developed by 
[Fleishman fc Kuznetsov (2010)^ which are used in practice as the forward fitting input. 
Then, having a fitting procedure resulting in fast and reliable finding of the true source 
parameters is exceedingly important. The problem here is that most of the minimization 
algorithms often find a local minimum of the normalized residual (or of the reduced chi- 
square), while the ultimate goal of the fitting is to identify the global minimum. So far. 
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Figure 3. Simulated image of the radio emitting loop source at 4 GHz as observed by an ideal 
radio heliograph with ^1" pixel size resolution (left). Example of the model (symbols) and fit (solid 
curve) spectra corresponding to one particular pixel of the image displayed in the left panel (middle). 
Model-to-fit comparison of the magnetic field (right). 

we determined that the simplex algorithm is very efficient in finding a local minimum. 
Then, it needs to be 'shaken' for the simplex solution to overcome any local minima and 
continue downhill towards the global minimum (a version of the stimulated annealing 
approach). Even when the algorithm performance is overall good, there is a non-zero 
probability that the algorithm fails to find the true solution in some pixels. We use 
post-processing to identify and flag/remove those pixels. 

3. Conclusions 

Modeling flare geometry, its full 3D visualization, and interactive adjustments to the 
user-specific needs are highly complicated tasks in themselves. Our modeling tools offer 
a united solution for these problems, which are widely applicable to various external 
data cube inputs and so offer a convenient framework for diverse studies of coronal 
magnetism, including flares and the active region magnetosphere. The modehng tools, 
computational libraries, and documentation are available via author's web page, see 
http://web.njit . edu/^gf leishm/. 

The outlined modeling efforts can only bring fundamental knowledge about flare/active 
region physics if used in conjunction with modern, high-resolution observations. Key ob- 
servations of the coronal plasma parameters can only be made by radio instruments that 
combine high sensitivity, temporal, spatial, and spectral resolution, which are unavailable 
now. A small part of the required science will be possible soon with the expanded OVSA 
instrument (anticipated operation of the upgraded instrument begins in fall, 2013) and 
the upgraded multi- wavelength SSRT. However, the full required capability has to wait 
until the fuh FASR (Gary 2003) has been built. 
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